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We use the three-body model and the sudden approxima¬ 
tion to compute angular correlations in high-energy fragmen¬ 
tation reactions of two-neutron halos on light targets. The 
contribution from one-neutron absorption by far dominates 
over that of neutron scattering. We use 6 He (n+n+a) and 
n Li (n+n-|- 9 Li) as examples and study the dependence of the 
predictions of this model on different physical assumptions 
and parameters. 

PACS number(s): 25.60.-t, 25.60.Gc, 21.45.+V 

Introduction. A new class of nuclear states, called ha¬ 
los, was discovered about ten years ago |l j. Halos are spa¬ 
tially extended states with small one or two-nucleon sep¬ 
aration energies and low orbital angular momenta. The 
systems can be understood essentially as few-body sys¬ 
tems where we can separate the coordinates into tightly 
bound intrinsic (core) and loosely bound external (halo) 
degrees of freedom. The characteristic properties, un¬ 
usual in the nuclear context, are related to the few halo 
degrees of freedom while the core can be assumed to be 
inert. Prominent examples are the two-neutron halos 
consisting of two neutrons and a core. The properties 
of halo systems have been remarkably well described by 
three-body models 01- The lesson is clearly that prop¬ 
erties related to the halo degrees of freedom (the core 
remains intact) can be deduced from three-body models. 

However, the recently measured neutron angular cor¬ 
relations in fragmentation reactions of Borromean two- 
neutron halos are apparently sensitive to both the reac¬ 
tion mechanism and the halo structure [Q . A satisfactory 
description is not available. Other previously measured 
neutron momentum distributions in high-energy breakup 
reactions have been reproduced in a three-body model 
where the sudden approximation and the final state inter¬ 
actions are decisive ingredients HQ] . Recently this model 
was extended to include both absorption and diffraction 
of the halo particles on the target and in addition also the 
absolute two-neutron removal cross sections were success¬ 
fully calculated ||. The purpose of this report is to ana¬ 
lyze the recently measured neutron angular correlations 
within this new model. We shall study the dependence on 
various physical parameters, compare to available data, 
see how far the model can go, possibly suggest new re¬ 
action mechanisms and predict yet unobserved neutron 
correlations. 

Model and method. The spatially extended three-body 


halo collides with a relatively small target at high en¬ 
ergy. Then the probability that more than one of the 
constituents interacts strongly with the target is very 
small. The differential cross section dcr is then to a good 
approximation a sum of three terms da W each describ¬ 
ing the independent contribution to the process from the 
interaction between the target and the halo particle i. 
This is the assumptions used in the classical formulation 
for a weakly bound projectile {tJ. We neglect the binding 
energy of the initial three-body bound state compared to 
the high energy of the beam. The reaction is then de¬ 
scribed as three particles independently interacting with 
the target as if each particle was free. 

The process is described as removal of one particle 
(participant) while the other two particles (spectators) 
both survive undisturbed. The participant is either ab¬ 
sorbed or elastically scattered by the target. The final 
state then consists of two independent subsystems, i.e. 
the two spectators and the target plus participant. The 
final state interaction between spectators is necessarily 
the same as in the initial three-body bound state which 
is described in agreement with available experimental in¬ 
formation [HQ . Since we only need to account for details 
of the scattered particle we employ a phenomenological 
optical potential where all other processes are included 
as absorption from the elastic channel ||. 

The coordinates, r^, R. and R', used to describe the 
reaction are sketched in Fig. Q. We denote the conjugate 
momenta by the corresponding p and use primes for the 
final states. We neglect the Coulomb interaction and 
assume that the target has zero spin. With the assump¬ 
tions of two independent subsystems in the final state, we 
must find both elastic and absorption halo-target differ¬ 
ential cross sections as products [Q of participant-target 
(elastic or absorption) cross sections and the averaged 
spectator overlap matrix element [Q, i.e. 

dP'dp' jk d p' 0z 

(PQi^Poi) 1 V- \M jm I 2 

dp' Qi 2J+1 Sjfc Si I 1 Sjk^2jk^2i I 

= 

dP ' d p’jk 

a ibs(P0i) 27+1 TjMs jk T. jk Y.i \ M s J 1 kE jk S i \ 2 > 


(1) 

( 2 ) 


where 


d 3 ^^ (PtH^Ppi) 

A4 


is the participant-target differen- 


1 








tial elastic cross section, <j^ is the participant-target 
absorption cross section, Ei and E' are spin projections 
of halo particle i before and after the reaction, Sjk and 
E jk = E j + Efc are total spin and projection of the halo 
particles j and k, M/ A £. e . is the overlap matrix el¬ 
ement between initial and final states of the spectator 
wave functions || and J and M are total angular momen¬ 
tum and projection on the beam direction of the initial 
halo wave function. 



FIG. 1. Sketch of the reaction and the coordinates used. 
The target is labelled by 0 and {i,j, k} label the particles 
within the three-body projectile. 


In the experiment the core and a neutron are detected 
in coincidence with velocities approximately equal to the 
beam velocity. Thus in agreement with the participant- 
spectators approximation one neutron reacts with the 
target without destroying or significantly affecting the 
motion of the core and the other neutron. The finite ex¬ 
tension of the projectile and the target therefore requires 
that only configurations where the participant is suffi¬ 
ciently far away from the spectators can contribute to 
the reaction. We account for this by omitting those ge¬ 
ometric configurations in the initial wave function where 
the participant ( n) is closer to the two spectators (n and 
c) than cutoff distances r nn and r nc which are treated as 
parameters. This shadowing effect substantially reduces 
the absolute values of the cross sections |||| . 

Observable and parameters. Recently a new observable 
was measured after fragmentation of a two-neutron halo 
nucleus on a carbon target. The projectile interacts with 
the target and the angular distribution of the relative 
momentum of the detected neutron-core system is mea¬ 
sured in a coordinate system with the 2 -axis along its 
center of mass momentum j(ij] . This process has a contri¬ 
bution of about 70% from absorption of the participant- 
neutron 0- The remaining about 30% arises from elas¬ 
tic scattering of the participant-neutron which leaves two 
neutrons in the final state. Subsequently the equal con¬ 
tributions (« 15%) from the spectator-neutron (absorp¬ 
tion like distribution) or the participant-neutron are mea¬ 
sured. The contribution from the interaction of the core 


with the target is expected to be negligibly small |5) . 

We shall first concentrate on the « 85% of the absorp¬ 
tion like distribution. Then the relative momentum is 
p j k and the center of mass momentum is the conjugate 
final state momentum p' ,- k of the coordinate r* jk, see 
Fig. [I]. The angular correlation is then computed from 
eq. d) by integration over all momentum coordinates ex¬ 
cept the angle 9 between these p' fc and p' - k . 

We consider the nuclei 6 He (n+n+a) and n Li 
(n+n+ 9 Li) with the wave function obtained by solving 
the Faddeev equations in coordinate space Q with the 
potentials from [[□]] . The resulting three-body wave func¬ 
tions have 88% of p 2 - and 12% of the neutron-core s 2 - 
configurations for 6 He and 20% of p 2 - and 80% of s 2 for 
11 Li. The binding energies and the root mean square 
radii are (0.95 MeV, 2.45 fm) and (0.305 MeV, 3.34 
fm), respectively. For the neutron-target interactions we 
use non-relativistic optical potentials || with the phe¬ 
nomenological neutron- 12 C parameterization EDAI-C12 
@ valid for a range of neutron energies from 29 to 1040 
MeV. We include 35 partial waves in the calculations. 

The binding energies and sizes of the initial states and 
the experimental neutron and core momentum distribu¬ 
tions are essentially reproduced with these parameters 
H||. Furthermore, shadowing parameters of 3 to 4 fm 
maintain or even improve the agreement of these re¬ 
sults for 6 He and in addition the available absolute two- 
neutron removal cross section is reproduced ||. 

Partial wave division. The initial three-body wave 
function contains in our examples both s 2 and p 2 rel¬ 
ative neutron-core configurations. Removal of one neu¬ 
tron leaves the remaining neutron-core system (specta¬ 
tors) in a mixture of s and p-waves, which after the ab¬ 
solute square and subsequent integration allow diagonal 
s and p-terms as well as an sp-interference term. The 
shapes of the resulting three angular distributions dif¬ 
fer substantially and the weighting, which of course is 
predicted by our model, is decisive. However, the exper¬ 
imentally preferred reaction mechanism might turn out 
to be different. To gain insight we show in Fig. |] the 
contributions from these different partial waves. 

The s-waves are angle-independent and insignificant 
for 6 He and dominating for 11 Li. The p-waves vary sym¬ 
metrically with angle and the asymmetric sp-interference 
terms change sign for cos 0 = 0 and would therefore not 
contribute to the total cross section obtained by integra¬ 
tion over the angle 6. Thus the angular correlation for 
6 He is essentially determined by the contributions from 
the p-wave and it is therefore symmetric. For n Li the 
angular variation is essentially due to the sp-interference 
and therefore very asymmetric. 

We also in Fig. |2| show the contributions arising from 
three different angular momentum projections along the 
neutron-core (spectators) center of mass momentum. For 
simplicity we assumed zero 9 Li-core spin (s c = 0) in these 
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computations. The difference between results for s c = 0 
and s c = 3/2 is visible but not substantial. The asym¬ 
metry is for both nuclei again due to the sp-interference 
appearing only for to/ = 0, which contains part of the 
p-wave and all the s-wave contributions. These terms 
with to/ = 0 exhibit a strong angular variation. The 
terms with to/ = ±1, arising entirely from p-waves, are 
symmetric due to the lack of sp-interference and vary 
relatively little with angle. 
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more visible, since the predominant removal of the dom¬ 
inating s-states now produce a more p-like structure, i.e. 
more symmetric but still varying with angle. 
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FIG. 3. The total absorption distribution for the same re¬ 
actions as in Fig. The results for different shadowing pa¬ 
rameters are compared with the available measured distri¬ 
butions |J. For comparison calculations without shadowing 
and without final state interactions are shown. All calculated 
curves are normalized to the same value. 


FIG. 2. The absorption contribution, computed in the pro¬ 
jectile center of mass, from the spectator-neutron (« 70%) to 
the neutron-core (spectators) angular correlation after frag¬ 
mentation of 6 He (left) and n Li (right) on 12 C at 240 MeV/u. 
The shadowing parameters are equal r nn = r„ a and 3 fm 
for 6 He and 4 fm for 11 Li. We show the contributions from 
the neutron-core (spectators) relative s and p-waves, their 
sp-interference and the three different projections on the 5 He 
or 10 Li center of mass momentum which is very close to the 
beam direction. To obtain the contributions from the individ¬ 
ual projections we assumed zero core spin of 9 Li (s c = 0) for 
11 Li. The total result for s c = 3/2 is also shown (right). The 
curves are shown with the relative normalization obtained 
from the initial wave function. 

Effects of shadowing. We compare in Fig. |] the mea¬ 
sured angular distribution with computations for various 
shadowing parameters. The first impression is that the 
computed and measured distributions are quite differ¬ 
ent. Both the asymmetry and the variation with angle 
are much larger for the computed curves. 

Shadowing is simply removal of unwanted geometric 
configurations in the initial wave function. We first only 
remove the part of the three-body wave function within 
a sphere around the core (long-dashed curves), which 
implies predominant removal of the relative neutron- 
core (spectators) s-states. We obtain an essentially un¬ 
changed distribution for (> He, since the s-wave only con¬ 
tributes marginally (« 12%) and all the p-waves are re¬ 
duced by the same amount. The result is a reduction 
factor independent of angle. For 11 Li this shadowing is 


If we instead only remove a sphere around the 
spectator-neutron (dot-dashed curves) we predominantly 
remove relative neutron-neutron s-states, which must be 
transformed to the neutron-core Jacobi coordinate sys¬ 
tem for the spectators. For 6 He the removed wave func¬ 
tion mostly consists of relative neutron-neutron s-waves 
with one nodes (hyperspherical quantum number K = 2). 
The transformed wave function then has a smaller rela¬ 
tive content of the dominating p-waves and the p and 
sp-contributions decrease relatively compared to the s- 
wave contribution with a smaller asymmetry as the con¬ 
sequence. For 11 Li both the removed and transformed 
wave function is mostly s-waves without node (K = 0). 
The reduction of the dominating s-wave contribution 
then produces a more symmetric p-like angular distri¬ 
bution. 

The final state interaction is essential to reproduce the 
observed narrow neutron momentum distribution Q. For 
the angular correlation the effect is significant but still 
substantially smaller than the difference to the measured 
distribution. The final state interaction preserves the 
total angular integrated contribution from each partial 
wave. The contribution from the neutron-core (specta¬ 
tors) relative s-wave is angle independent. However, in¬ 
cluding the final state interaction and combining with 
the p-wave contributions the asymmetry for 6 He or n Li 
increase or decrease, respectively. 

Dependence on spectator excitation energy. We have 
assumed that the reaction is dominated by one-neutron 
absorption while the remaining neutron-core (spectators) 
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system is undisturbed. The argument is that this high- 
energy process must be very fast compared to the time 
scale of the intrinsic motion of the halo particles. The 
reaction has occurred before the remaining halo particles 
can change their relative motion and for example select 
an excited state which in the present cases must be two- 
body resonances or virtual states. This is the sudden 
approximation used so successfully for many other ob¬ 
servables pppl- In this picture all continuum states 
are populated with the probability of occurrence in the 
initial wave function. However, it is conceivable that a 
different continuum population arises in these reactions 
and in particular when the beam energy is substantially 
reduced below the values in the present experiments. 
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FIG. 4. The total absorption distribution for the parame¬ 
ters and reactions in Fig. H. The results are given for vari¬ 
ous relative neutron-core (spectators) excitation energies E na 
and E nc = p(f c /2/x. All calculated curves are normalized to 
the same value. 


To investigate the sensitivity of the angular correla¬ 
tion to the continuum population we show in Fig. [| the 
results when only a state with a given relative neutron- 
core (spectators) excitation energy is populated in the 
final state. This energy is varied from zero and across 
the energies of the lowest resonances Jll||. We find for 
6 He that all the distributions vary too much with angle. 
The lowest energy of 0.2 MeV exhibits a large asymme¬ 
try indicating a relatively large admixture of s-waves still 
far from the measured distribution. The higher energies 
all produce similar distributions even when the energy 
matches the p-resonance at 0.77 MeV. The continuum 
two-body states are essentially p-waves due to the low- 
lying resonance and the distributions in fact resemble the 
p-wave distribution from Fig. |2j. Thus, no weighted av¬ 
erage of these continuum states can reproduce the mea¬ 
sured distribution. 

For 11 Li the low-lying virtual s-states influence the an¬ 
gular distribution. We find a rather flat and asymmetric 
distribution at low energy, where the s-waves dominate 


completely. As the excitation energy increases the p- 
waves contribute more and more. This produces an in¬ 
creased asymmetry due to the interference term and at 
higher energies the p-wave contribution is clearly pulling 
the distribution towards symmetry. 

Contribution from the participant-neutron. If the 
participant-neutron is scattered by the target instead of 
being absorbed it may still be detected and contributes 
then about 15% to the measured cross section. The es¬ 
timate assumes that this neutron arrives within the for¬ 
ward angle where the detection takes place. This rela¬ 
tively small contribution is for technical reasons only esti¬ 
mated approximately. We first approximate the motion 
of the center of mass of the neutron-core (participant- 
spectator) after the reaction to be in the direction of the 
beam. Then the participant-neutron and spectator-core 
relative momentum is approximated, as for an infinitely 
heavy core, by the momentum of this participant-neutron 
relative to the projectile center of mass. The heavier the 
core and the higher the beam energy the better the ap¬ 
proximations. Finally we obtain the contribution from 
eq.®) by integration over all momenta except the angle 
between these two momenta. The result is shown Fig. [^. 



COS(9) COS(0) 

FIG. 5. The angular distribution for the parameters and 
reactions in Fig. pj. The insets show the two contributions 
from the scattering process with two neutrons and a core 
in the final state. Each of these contributions amounts to 
about 15%, i.e. when the participant (dashed) and the spec¬ 
tator-neutron (solid), respectively are measured. The solid 
curves in the external parts give the total contributions when 
the approximately 15% from the participant-neutron is com¬ 
pletely neglected, i.e. we show the sum of the solid curve (« 
15%) from the insets and the total contribution (« 70%) from 
Fig. The dashed curves in the external parts show the sum 
of the solid curves (« 85%) and the dashed curves (« 15%) 
from the inset. 

The contribution from the participant-neutron is al¬ 
most symmetric and peaked at cos 0 = 0, see the inset of 
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Fig. |. This behavior is in striking contrast to the contri¬ 
bution from the spectator-neutron which has a minimum 
at cos 9 = 0 for 1 2 3 4 * 6 7 8 9 He and is rather flat for 11 Li. The peak 
corresponds to a preferred direction perpendicular to the 
beam direction as obtained for forward scattering of a 
neutron on a target expressed in the rest system of the 
neutron. Such a distribution is broadened by the intrin¬ 
sic motion of the scattered participant-neutron with re¬ 
spect to the center of mass of the two spectator-particles. 
Thus, the higher the beam energy and the smaller the 
halo binding energy, the narrower the distribution. 

For 6 He the contribution from the participant-neutron 
only changed the total distribution very little, see Fig. ||. 
For 11 Li a peak appears due to the rapid angular varia¬ 
tion of the contribution from the scattered participant- 
neutron compared to the almost constant background 
from the dominating s-waves, see Fig. ||. Before a di¬ 
rect comparison with measurements it is worth keeping in 
mind that this result is obtained first as a crude estimate 
and furthermore with all participant-neutron contribu¬ 
tions included independent of possible additional exper¬ 
imental selection. 

Conclusion. We computed the recently measured an¬ 
gular correlation for breakup reactions of halo nuclei 6 He 
and 11 Li, i.e. the halo nucleus interacts at high energy 
with the target and the angular distribution of the rel¬ 
ative momentum of the detected 5 He and 10 Li neutron- 
core system is measured in a coordinate system with the 
z-axis along the center of mass momentum. We used a 
model which successfully describes essentially all other 
three-body observables for such systems. 

We investigated the contributions from different par¬ 
tial waves of the relative motion of the detected neutron- 
core system and found constant, symmetric and asym¬ 
metric distributions arising respectively from diagonal s, 
diagonal p and sp-interference terms. The zero angular 
momentum projection on the direction of the 5 He and 
10 Li momenta produces by far the largest variation with 
angle. 

Exciting continuum neutron-core states of a given en¬ 
ergy in the reaction with the subsequent decay of these 
states might be a possible reaction mechanism which 
however is not supported by the present analysis. The 
energy is too high to allow time to select specifically for 
example the two-body resonance states and agreement 
with the available data is furthermore not improved. 

Two processes are possible, i.e. absorption or scatter¬ 
ing of one neutron (participant) while the other neutron 
and the core (spectators) continue undisturbed. The de¬ 
tected neutron-core system may consist of either partici¬ 
pant or spectator-neutron and the resulting angular dis¬ 
tributions are qualitatively different. The process where 
the participant-neutron is detected contributes by less 
than the 15% predicted by the optical model, since only 
a narrow cone around the beam direction is selected in 
the experiment. 


The correlations for 6 He and 11 Li are qualitatively dif¬ 
ferent and the computed distribution vary too strongly 
compared to the data for 6 He. Provided the experimental 
correlation is correct this indicates that a subtle reaction 
mechanism is at work. From the present investigation 
this can be either a missing constant background for ex¬ 
ample due to suppression of p -waves or a preferred selec¬ 
tion of s-waves or a relative suppression of zero angular 
momentum projections. 
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